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I. INTRODUCTION 



Experimental studies of hard exclusive processes, in particular, the production of dijets 
at midrapidities, has been recently performed at Tevatron [l[ and will be done at the LHC 
in the near future @, 0]. 

It is expected that the contribution of the gluon pairs to the exclusive hard dijets produc- 
tion dominates over that from quark-antiquark pairs. Martin, Ryskin and Khoze proposed 
a QCD mechanism of exclusive digluon production 0| . In certain regions of the phase space, 
the process pp — > p(gg)p (similarly to qq production) is dominated by the non-perturbative 
region of gluon transverse momenta, and even perturbative ingredients like the Sudakov form 
factor are not under full theoretical control [6] . The problem becomes even more pronounced 
when considering the irreducible backgrounds in central exclusive production (CEP) of Higgs 
boson originating from the direct exclusive bb pair production in a fusion of two off-shell 
gluons. In particular, in Ref. 0, 0] it was shown that the central exclusive production of bb 
jets at the LHC, may noticeably shadow the corresponding signal of the Higgs boson in the 
bb decay channel. Along with unknown NLO corrections to the fc^-dependent hard subpro- 
cess amplitude g*g* — > jj (in particular, the NLO contribution from the rescattering of the 
final state gluons into qq pairs can be potentially important), this may lead to problems in 
experimental identification of the Higgs boson. 

Recently, it was shown in Ref. [§] that the first LHC measurements of the exclusive dijets 
would significantly reduce the theoretical uncertainty for the central exclusive Higgs boson 
production. This makes the process under consideration especially important from both 
theoretical and experimental points of view. 

Such a process has been recently investigated in detail in Ref. 0], and fairly good descrip- 
tion of the Tevatron data has been achieved. We would like to extend such an analysis, both 
analytically and numerically, by analyzing separate contributions from different final gluon 
polarisations, various kinematical regions of the 4-particle phase space, which are important 
for future LHC measurements, and the theoretical uncertainties related with different choice 
of UGDFs and the factorisation scale. Similarly to Ref. we shall limit ourselves to the 
lowest-order QCD calculation, and postpone the analysis of the higher order contributions 
for a separate study. 

Compared to the previous studies, we would like to perform an estimation of the pro- 
cess when one gluon is emitted from a one t-channel gluon of the QCD ladder and the 
second gluon is emitted from the other t-channel gluon. This contribution was discussed in 
Ref. p] as potentially sizable. In this work, we present the first numerical calculation of this 
contribution. 

Recently, the calculation of the exclusive production of quark-antiquark dijets has been 
performed for heavy cc T(| and bb 0, 0] pairs. Here, we extend our previous analysis and 
present the calculation including both contributions from light (w, d, s) and heavy (c, b) 
quark/ ant iquark jets compared to the one from the gluonic jets. 

A high precision measurement of exclusive bb pair production is required for central 
exclusive Higgs production measurements 0, 0]. It is, therefore, instructive to estimate the 
reducible background for Higgs CEP at LHC coming from misidentification of gluon jets. 

Another interesting point, which we would like to investigate here, is the role of the 
gluon reggeization in the exclusive gluonic dijets production. For this purpose, we employ 
the formalism of the quasi-multi-Regge kinematics (QMRK) with the Lipatov's nonlocal 
vertices for the triple-gluon coupling [ll|, and perform a numerical comparison with the 
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standard pQCD calculation (with standard gluons) [6|. 

This paper is organized as follows. In the second section, we present the standard exclu- 
sive diffractive amplitude when both gluons are produced from the same t-channel gluonic 
line as well as the amplitude when one of the gluons is emitted from one line and the 
second gluon from the second line. In the third section, we briefly remind formulae for 
quark-antiquark dijets with arbitrary quark mass. The fourth section contains discussion of 
unintegrated gluon densities and model assumptions. In the Results section, we present pre- 
dictions for various differential distributions and compare our results with the available CDF 
data as well as discuss corresponding theoretical uncertainties. We also present predictions 
for future studies at LHC Finally, the summary and conclusions close our paper. 



II. DIFFRACTIVE AMPLITUDE OF THE EXCLUSIVE GLUON PAIR PRO- 
DUCTION 



In this analysis, we apply the QCD mechanism for the central exclusive production, 
proposed by the Durham group (referred to below as the KMR approach) in Ref. [l2j . In 
Fig. [TJwe show typical contributions to the exclusive gluon pair production. In the important 
limit of high-p t jets, the diagram B is suppressed by an extra hard propagator. Such a limit 
has been considered in detail in Ref. [f|. However, at relatively small gluon p^'s, the diagram 
B may become sizeable and a reliable numerical estimation of its contribution is required. 
In this paper, for generality, we would like to calculate both contributions in all potentially 
interesting regions of the 4-particle phase space. 
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FIG. 1: Typical diagrams for the exclusive gluon pair production in exclusive double diffractive 
pp scattering through the gluon- gluon fusion subprocess g*g* — > gg (A) and the 4-gluon fusion 
subprocess g*g*g*g* — > gg (B). 



Momenta of intermediate and final state gluons are given by the following Sudakov de- 
compositions in terms of incoming protons momenta 

?i = X\Pi + Qi±, q2 = x 2 p 2 + q 2 ±, Qo = x'iPi +x' 2 p 2 + <?o_l ^ Qo±, x'i,2 < £1,2(2.1) 

P3 = 0iPi + am + h±, p A = 2 pi + a 2 p 2 + k 2± . (2.2) 
In forward scattering limit, we then have 

ti, 2 = (pi,2 ~ p'i, 2 ) 2 = p'\/2± -> °> 9o_l ^ ~qi± ^ 42_L • (2.3) 



P i/2± -> °> 9o± - -qi± - Q2± ■ 

The Mandelstam invariants in the two-gluon fusion in the limit of high-p^ gluon jets |k| 
|p 3 | ~ |p 4 | ^> |q | can be written as 0] 



M 



99 



°99 — 



; (/3l+/3 2 ) S 
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(2.4) 
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Here and below, we use notations for the transverse 2-momenta in bold face style. 

Let us first consider the explicit derivation of the diffractive amplitude shown in Fig. 
as an example. Starting at the parton-level process and applying the cutting rules the 
imaginary part of the one-loop partonic amplitude of the gluon pair production with a fixed 
color indices hi and b 2 can be calculated in the forward limit as (for similar derivation of 
Higgs CEP amplitude, see Ref. 0) 

2 2 (2tt) ; 

x / 4^2 d ^ dx 'MPi - qo) 2 m P 2 + go) 2 ) ■ Win, ri)V^(n, -q 2 ) ■ e^AiX (A 2 ) , 

where the first factor 1/2 comes from the cutting rule, factor 2 comes due to two identical 
contributing diagrams (emission of jets from the first and second t-channel gluon line), factor 
1(2^)2 comes from the phase space in the loop integration, P 12 are the effective Reggeon- 
Reggeon-gluon (RRG) verticies in the quasi- multi-Regge kinematics (QMRK), corresponding 



lmM parton = _ . 8 _ T^Tf^-fif^f*** • {2 9 s )\ Pl p 2 )p lp p 2(J • {i 9s f (2.5) 



to the kinematical configuration with /3i 3> /3 2 , o>i <C a 2 pjj. Here, we apply the eikonal 
approximation for the quark-gluon verticies in the proton defined by 2g s T^pi^\\'i where r a 
are the Gell-Mann matricies, g s is the QCD coupling, and 5w appears due to the fact that 
soft gluons cannot change quark helicity in an energetic proton. 

For the color singlet production (since in- and outgoing protons are in the color singlet 
state) there is no color transfer from quark lines between protons, i.e. color indices in initial 
and final quarks are the same, i.e. i = k and j = I. Then, the color averaging in each quark 
line leads to the substitution 

a a ci c 2 ^ 

irn'jn'mk nl 4_/\T 2 

in the diffractive amplitude (12. 5ft . The appearance of 6 ClC2 here automatically gives rise to 
the projection of the produced gg-pair onto the color singlet state. 

Further, in order to go over to the hadron level, one has to absorb the factor 

C F a s = N 2 C - 1 a s 

into a definition of the unintegrated (q±- dependent) gluon distribution function (UGDF) 
along each proton line as required by the underlying fc^-factorisation approach [13| . A reli- 
able model for generalized off-diagonal UGDFs used in the current analysis will be discussed 
in some detail below. 



Due to the gauge invariance and the factorisation property of the hard RRG vertices [11 
we get 



Pi^r (-)^r (-) = s -rn + P n-v^{...)vi^{...) -+ ^i(-)^(-) > ri = qi _ P3 

4 4 Y\ 

where n 1 * 1 = pi 2 /Ep ms , Ep ms = y/s/2, and Cf 2 are the nonlocal RRG couplings defined as 

H 

V2) = ( A - m) _ f2 ( ai _ g _ (V1± + V2±r , 

C^( Vl , v 2 ) = pi (fr - ^l) - j4 (a 2 -^f)- (vxx + v 2± r ■ (2-7) 
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Finally, the contributions of diagrams Fig. [T] (A) and (B) to the diffractive amplitude 
J\A" = M. A + M. B for the central exclusive gg (with external color indices a and b) dijet 
production pp — > p(gg)p read 



A4£(Ai,A 2 ) 



is A 



>ab 



d 2 q.o 



ff(qo,qi)ff(qo,q2)-e;(X 1 )et(X, 



X 



C^qurJC^-q,) Cffa, r 2 )Q(r 2 , -g 2 ) 

o T o 



A<£(Ai,A 2 



-is .4 



/f(«i,«3)/f(«2,«4)-6:(A 1 )e*(A 2 ) 



Cf (/C 1( -K 2 )C 2 (k 3 , -K 4 



Av-^ 



X 



(2.8) 



(2.9) 



where ^4 = 2n g s /Cp, the minus sign in .M comes from the difference in colour factors, 
/° ff (t>i, t> 2 ) is the off-diagonal UGDF, which is dependent on longitudinal and transverse 
components of both gluons with 4-momenta v\ and t> 2 , emitted from a single proton line, 
and 

Then the matrix element squared for the exclusive diffractive production cross section 
can be written in the standard way 



\M\ 2 = Y,M ab (M ab ) 



(2.10) 



a,b 



summing up over all possible color singlet combinations of the final gluons. 

The integration in the case of the emission from different t-channel gluon lines (diagram 
B in Fig. [T|) can be made symmetric with respect to both protons by the following equivalent 
transformation of the integral measure 



d 2 K,i 



d K-d k+5 (k. 



Pa) 



(2.11) 



where k_ = K\ — k 2 and n + = K\ + k 2 have been introduced. In practice, such a transfor- 
mation is convenient in numerical calculations below. 

When the p_i_'s of the final jets are sufficiently large, the contribution of the diagram B 
should vanish much faster than that of the diagram A due to an extra propagator suppression 
in the amplitude (see Eq. 12. 9f) . Moreover, the hard transverse momentum flow through the 
proton remnant would disturb it too much such that it becomes less likely to combine it back 
to an exclusive proton-like state after hadronisation, which should be reflected in an extra 
suppression by the UGDFs behavior at large gluon p±. Such arguments lead to a conclusion 
that the diagram B can be sizeable only at relatively small jet p_i_'s, but the large invariant 
mass of the gg dijet system (i.e. at large rapidity difference Ay = \yi — y 2 \ between two jets 
at the edges of the central detector). Numerical estimation of such a contribution could be, 
therefore, important when the statistics on diffractive dijets production at LHC becomes 
sufficiently large. Potentially, it could even be singled out and its features could be tested if 
there is a region in the phase space where it may dominate. 

For the emission of both gluons from the same t-channel gluon line (the standard CEP 
process given by diagram A in Fig. [[]) we have typically: q±± ~ g 2 j_ p%± ~ p^±. The 
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integration over screening gluon transverse momentum is limited to rather small qo± ~ qi/2± 
(in the forward limit). 

The kinematical situation for the diagram B is different. Typically, in this case either k%± 
is large and of the order of p%± and k 2 ± is small, or vice versa — K\± is small and k 2 ± is large 
and of the order of p^±. The integration over k_ extends to large values, which means that 
typical transverse momenta of gluons in the impact factors are large. So, technically, when 
using grids for UGDFs we have to do it separately for both situations. Also, the kinematical 
structure of UGDFs are very different in diagrams A and B. This issue will be discussed in 
detail below. 

In Eqs. (12 .8p and (I2.9p . e*(Ai) and e*(A 2 ) are the polarisation vectors of the final state 
gluons with helicities Ai, A2 and momenta ps, p±, respectively. They can be defined in the 
gg rest frame with z axis along the proton beam as 

e*(Ai) = ^(0, Ai cos 9 cos ip — isimp, Ai cos^sm-?/? + icos^, — A x sin^) , 

e*(A 2 ) = ^(0, — A2 cos 9 cos ip — i simp, — A 2 cos 9 simp + % cosip, A2 sin 9) , (2.12) 

a/2 

such that e M (Ai)e*(A 2 ) = -<5 Al ~ A2 and e*(Ai)p3 = e*(A 2 )p4 = 0. In this frame, momenta of 
protons and final-state gluons are 

Pi = t|(1> 0, 0, 1), p£ = %1, 0, 0, -1), 

P3 = E g (l, sm9cosip, sin 9 sin ip, cos9), (2-13) 
p\ = E g (l, —sm9cosip, — sin 9 siiaip, —cos9), 

so that the proton and gluon energies Ei )2 , E g and the polar angle of a gluon jet 9 w.r.t. 
the z-axis are defined as 



r \/Sgg E x E 2 , 0\ — 2 . a 2a//3i/3 2 

Eg = JL —— = — (pi + p 2 ) = — (ofi + a 2 ), cos9 = - — r - B -, sm9 



a/2 ' y/T ' 01+02 01+02 

In the high-p^ limit and at central rapidities of jets, the gg dijet rest frame, introduced 
above, becomes identical to the initial protons c.m.s. frame, which we use in actual numerical 
calculations below. Then the diffractive amplitude f !2.8j) reduces to the standard expression 
with gg — > gg hard scattering amplitude initially derived in Ref. j6[ 

Sab f , 2 /f (g O,gl)/f ( g0,ga) 
q q l q 2 

(-2AtA;)| qi ||q 2 |e- iA ^ 1+iA ^ 2 • A(\l\* 2 ->• AiA 2 ) , (2.14) 

where 0i j2 and <fi are the azimuthal angles of the fusing gluons gi j2 and of the dijet production 
plane, respectively, \\, A 2 are the helicities of fusing gluons with momenta q±, q 2 , respectively, 
and the nonzeroth helicity amplitudes 



™ ab - 2tA^ J q q2q2q2 ^.e 



A(++ ++) = A( ► — ) = 1 

A(+- -> +-)=A(-+ -> -+) 



A(+ ► -+) = A(-+ -> +- 



u 2 

U 99 
t 2 

s 2 ' 

4 9 '9 
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with Mandelstam invariants denned in Eq. (12. 4p . from which we see that in the QMRK limit 
amplitudes A^XIX^ — > A1A2) reduce simply to <5a*A!<5a*a 2 - 111 the forward limit, provided by 
Eq. (12. 3p . we have 02 — 4>i + n, so the integral in Eq. (I2.14p 




-i(Ai— Az)^i 



survives only when Ai = A2, i.e. when the gg dijet is produced in the J z = state, which 
corresponds to the well-known J z = selection rule in the central exclusive production 
processes [151 ] . We will discuss the subleading corrections to this rule below when presenting 
the numerical results for different gluon polarizations. 

Below, for consistency, in order to see effects of the gluon reggeization and subleading 
corrections to the high-p^ limit we numerically compare both versions of the diffractive 
amplitudes - in the QMRK approximation (I2.8P and in the standard approach (I2.14p . as 
well as estimate subleading contribution to the observable signal from the amplitude B with 
symmetric gluon couplings (12. 9p . The bare amplitudes above are subjected to absorption 
corrections which depend on collision energy and typical proton transverse momenta. We 
shall discuss this issue shortly when presenting our results. 



III. QUARK- ANTIQUARK DUETS PRODUCTION 

Let us consider now the contribution of the quark/ ant iquark pairs to the observable signal 
of the exclusive dijets production. The hard subprocess amplitude for the qq pairpro duct ion 
via off-shell gluon-gluon fusion was previously discussed in detail in Refs. [71, [sj] . Here we 
would like to list the relevant formulae only, which will be used in numerical calculations 
below. 

The amplitude of the exclusive diffractive qq pair production pp — > p(qq)p reads jl] 

M„iX u X,) = is-^jjp^, f^V?T / f(g °'^£ {q °' g2) , (3.1) 



Here, the vertex factor V^ 2 is the production amplitude of a pair of massive quark q and 
antiquark q with helicities X q , Xq and momenta ps, p 4 , respectively. It is given by the 
following general expression 

M qq±VN^ \ ?r±-2(p 3 gi) gr ± -2(p 4 gi) / 

(3-2) 

In analogy with Eq. (12. 2p . one can introduce the Sudakov expansions for quark momenta 



as 



P3 = Avi + AV2 + P3±, Pi = Avr + X \V2 + Pi± (3.3) 

leading to 

x i,2 = A,2 + s i,2> x i,2 = f ^^e ±v \ x\ 2 = r ^^e ±v \ m 2 l/2± = m 2 q + |k 1/2 _L| 2 , (3.4) 
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in terms of quark/antiquark rapidities y±, yi and transverse masses mu_, m 2 j_. The only 
difference of the quark/antiquark fractions x\\ from that of the gluons is that they are 
dependent on the quark mass m q . 

It is convenient to fix the c.m.s. frame of the qq pair with z axis along the proton beam, 
so p 3 = — p 4 = k and P34 = M qq /2. The gluon and quark transverse momenta (with respect 
to the proton beam) in the polar coordinates are then defined as 

qo_i_ = qj_(cosij;, simp), p 3 j_ = — p4j_ = A;_i_(cos x, sin x), 

respectively, and 



k\=E, 



^(xi + xiy-ixi-xiy 



1 1 

k = E Xl ~ Xl 



A;2 + fc2 = Eql . (3.5) 



In these notations, the helicity amplitudes g*g* — > q(Xi)q(\2) can be written as (for more 
detail, see Ref. [sj]) 

2|k|gj_ ^|k| cos(^ — x) — ik z sin(ip — x) j + M qq kj_ (k z cos(2?/> — 2x) — 



|k| sin(2'0 — 2x^ 



/ 



M U k l + 9± + + 4M qq -k ± q ± k z cos(^ - x) 



2 ^igi(l + cos(2^ - 2x)) + ^ 



(3.6) 



V++ = -2C e" 



£r[cos(2^ - 2x) + |k| 



/ 



MUkl + Q 2 ± + ml) + 



iM qq k±q±k z cos(^ - x) - 2k\q]_{\ + cos(2^ - 2x)) + q\ 



(3.7) 



where m g is the quark mass, and the normalisation factor is C = 2g 2 S ClC2 / 'yfW c . Below, we 
will use the expressions (13.61) and (13 .7p in calculations of the quark jets contribution to the 
central exclusive dijets production at Tevatron and LHC energies of both light (u, d, s) and 
heavy (c, b) quarks. 



IV. OFF-DIAGONAL UNINTEGRATED GLUON DISTRIBUTIONS 
A. Emission from the same i-channel gluon 

The off-diagonal unintegrated gluon distribution in Eq. (I2.8p . where longitudinal momen- 
tum fractions satisfy the strong inequality x' <C x, is calculated in the forward scattering 
limit qo± ~ qi/2± according to the Kimber-Martin-Ryskin (KMR) prescription [l6|, 17 



fif 2g (xi,2,x',q 2 /2± ,ql ± ^ 2 ;t) = R g( x ') j^2 ( x lM x i,2, ql)S(ql, /i 2 )) \ q l= q f /2X -F(t) , (4.1) 



which was used e.g. in exclusive dijets studies in Refs. |9[ and leads to a reasonable de- 
scription of the Tevatron data [lj. In the equation above, X\p and ? 2 / 2 j_ are longitudinal 
momentum fractions with respect to the parent proton and transverse momenta squared of 
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the active gluons q x ± and q 2 ±, respectively, x' ~ qoj_/y/s and q% L are the same variables 
for the screening gluon q ±. If we assume that at small x: xg(x) = N g x Xa the skewedness 
parameter can be expressed in terms of the X g as (l8| : 

2^ r(A J + 5/2) 

R ° = ^ r(A 9 + 4) • (4 ' 2) 

We will take i? 9 ~ 1.2 in practical calculations. The function ^(a;, g^_) in Eq. (14. ip is the 
collinear DGLAP gluon distribution, S(q\,fi 2 ) is the so-called Sudakov form factor and 
the nucleon form factor in the forward limit F(t) = exp(6t/2) with the slope parameter 



b ~ 4 GeV 2 [19j describes the coupling of the gluonic ladders to one of the nucleon lines 
(see Fig. d) . 

The range of the integration over <i 2 qo is formally limited by the existence of the DGLAP 
gluon PDF xg(x,q"^). In explicit calculations we use the next-to-leading order CTEQ6 
collinear distributions [20] for which xg(x,fi 2 ) parametrization works well down to quite 
small /2q ~ 0.4 GeV 2 . This way we cut off the region below this starting scale /i 2 ,. This 
has not practical consequences for production of relatively large invariant masses as these 
regions are also numerically suppressed by the behavior of the Sudakov form factor for small 
q\. In order to exhibite uncertainties of our numerical results related to the collinear PDFs 



we also use GRV94 |2lJ, GJR08 |22j and MSTW08 |23| distributions. 

The longitudinal momentum fractions of the fusing gluons entering Eq. (14. ip are calcu- 
lated as 

P3_|_ P4_|_ 

x x = —= exp(+y 3 ) + — exp(+2/ 4 ) , 
V s V s 

x 2 = ^exp(-y 3 ) + ^exp(-y 4 ) (4.3) 
V s V s 

from the transverse momenta P3/4,j_ and rapidities 3/3,4 of the final gluonic jets. 



B. Emission from both t-channel gluons 

The kinematics of the diagram B in Fig. [T] is different from that for the diagram 
A. Here, the off-diagonal unintegrated gluon distributions f° G (xi, X3, n\ , «§, /1 2 , and 



, x 4 , k|, k 2 , /i 2 , H2) should be evaluated at X\ ~ x 2 and 23 ~ X4. In general, such 
objects are not well known and were not discussed so far in the literature. We calculate the 
longitudinal momentum fractions of the fusing gluons in the considered kinematical domain 
as follows 

P3.L , . v P4± ( v 

x x ~ — exp(+y 3 ) , x 2 ~ —= exp{-y 3 ) , 



S y/S 

x 3 ~ exp(+y 4 ) , x 4 ~ exp(-y 4 ) . (4.4) 
V s V s 



As a first approximation, one could try to use a symmetric factorized prescription for the 
off-diagonal UGDFs, which was successfully used before for the exclusive production of \c 



mesons in Ref. 24 



f° (xuxa^^n^fi^^t) = J/ g (xi,K 2 ,/i 2 )/ 9 (x 3 ,K 2 ,/i 2 ) ■ F(t x ) 



/° (x 2 ,x 4 ,K 2 ,K 4 ,ii 1 ,n 2 ;t) = J fg(x 2 ,Kl,fij)f g (x 4 ,Kl,fil) ■ F{t 2 ) . (4.5) 
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Above unintegrated diagonal distributions include Sudakov form factors in the same way as 
in the KMR UGDF (14. ip . Since for the jet production in the diagram B pz± > Ki±, k 2 ± 
and p4j_ > k>3±, K4j_ in most cases, so a physically reasonable choice of scales in the scale- 
dependent UGDFs would be /ii = p 3 j_ and /z 2 = p 4 j_ or ^ — // 2 = Mjj. We adopt these 
simplest choices since we do not know the exact evolution of the Sudakov form factor in the 
considered kinematical domain (see, also a discussion of this issue in Ref. \$\). In this case, 
when the gluon q± becomes bigger than the scale /x, we take simply S(q\,fi 2 ) = 1. 

Our prescription in Eq. (14. 5 p does not use the fact that in the considered process both 
gluons are outgoing (emitted) . In the collinear approach this corresponds to the ERBL 
kinematical region [25j] where 

ii , H x 1 + x 2 Xi-X 2 

In our case of central dijet production typically both x\ and x 2 are small that is also x 
and £ are small. In this region the collinear off-diagonal distributions H(x, £, fi 2 , t) can be 
estimated in a model independent way (l8| . 

The discussion above suggests therefore another prescription for UGDFs in this special 
kinematical case: 



ff(x, x', k 2 tl k' 2 t , /i 2 , t) = R coll (x, x'- /i 2 , t = 0) • yjf g (x,kl^)f g (x,k'l^) ■ F(t) , (4.6) 

where x = p? = /i 2 ~ /i 2 and f g are standard diagonal unintegrated distributions as in 
e.g. Ref. 16J, ll7| . Here R co ii is the ratio of collinear off-diagonal distributions in ERBL to 



DGLAP region: 

H z RBL ( x ^;fi 2 ,t = 0) 
Reafaxw ,t = 0) = H nc L AP {x ^.^ t = 0) ■ ( 4 - 7 ) 

Assuming that at small x: xg(x) = N g x~ Xg in the limit of small x and £ the off-diagonal 
distribution H g (x,£,t) can be expressed in terms of X g as: 

For our estimates here X g is a crucial parameter which is not completely well known. In the 
double logarithm approximation at small values of x: 



W^ log log @)' (4 - 9) 

However, the gluon distribution at x < 10~ 4 and small factorization scales is poorly known 
(see e.g. a discussion in Q). In consequence applicability of the double logarithmic formula 
( 14. 9 p is not obvious and not well justified. Therefore we will treat X g as a free parameter. 
In general, it can be dependent on the scale of the problem (transverse momentum of the 
jet). To demonstrate uncertainties we shall show results for \ g = 0.2, 0.4, 0.6 and 0.8. 

In Fig. |5]we show the ratio R co ii for different values of X g . The ratio strongly depends on 
the value. We observe a strong enhancement on the diagonal. The larger X g the stronger the 
ratio. The ratio quickly drops off diagonal. This have consequences for rapidity distributions 
of jets, in particular their correlations, as is discussed in the next section. 
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FIG. 2: The ratio from Eq. fl4.7j) as a function of xi and x 2 for different X g = 0.2, 0.4, 0.6, 0.8. 



V. RESULTS 



Let us start presentation of our results. Having the amplitude M. for the pp — > p(jj)p 
where j = g,q(q), defined in Eqs. (12. 8p . (12. 9 p and (13. ip . we can calculate the corresponding 
2—7-4 cross section as 



f d 3 p 1 d?p 2 d 3 p 3 d 3 p A , 4 (4) , . 



a pp^pjjp 



In what follows, we adopt a convenient choice of the phase space variables of the integration 



relevant for exclusive diffractive processes elaborated in Ref. [26| . 

Before we go to the description of experimental data and presentation of all contributions 
let us concentrate for a while on the contribution of the diagram B mechanism of hard 
digluon production. In Fig. [3] we show distributions of the gluonic jets in pseudorapidity 
and transverse momentum of the jet. The results have been performed for different values of 
X g using formula (14. 6 j) with R co u based on H g from formula (14.81) . The rapidity distribution 
strongly depends on the value of X g . The dependence is stronger for smaller transverse 
momenta, i.e. smaller x's. 

In FigJUwe show similar distributions for the dijet invariant mass (left panel). There is 
a stronger dependence on X g at small invariant masses. In the region of the Higgs boson 
invariant mass of Mjj = 120 GeV there is a factor three uncertainties of the cross section. 
For comparison we show result (dashed line) obtained from a naive prescription (I4.5p . When 
the naive prescription gives large cross sections at large invariant masses the improved 
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FIG. 3: Distribution in pseudorapidity (left panel) and transverse momentum (right panel) of 
gluonic jets for the contribution of diagram B. We present results for different X g = 0.2, 0.4, 0.6, 
0.8. The results strongly depend on the value of X g . For comparison we show result (dashed line) 
obtained from the naive prescription f|4.5j> . 



prediction drops quickly with invariant mass. For completeness in the right panel we show 
disitrbutions in rapidity difference between jets. While the distribution obtained with naive 
prescription extends up to large rjdiff, the distribution obtained with improved calculations 
is concentrated at small values of rjdiff. 




FIG. 4: Distribution in dijet invariant mass (left panel) and in pseudorapidity difference between 
jets (right panel) for the contribution of diagram B. We present results for different X g = 0.2, 0.4, 
0.6, 0.8. The results strongly depend on the value of X g . For comparison we show result (dashed 
line) obtained from the naive prescription (|4.5p . 

Now we shall include and discuss all contributions. First we shall discuss numerical results 
obtained at the Tevatron energy and we wish to compare our results to the existing CDF 
collaboration data [1]. This means that in the following we include the CDF experimental 
cuts. In Fig. we show the integrated cross section as a function of lower cut on E T . 
Following Ref. [6[, we assume the relation between Et and the jet transverse momentum p± 
as Et = 0.8 p±. This approximate relation could be checked in the future by performing 
full simulation of jets including hadronization. We show results for the digluon diagram 
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A (solid curves), as well as for the diagram B (left panel, long dashed curve) and for the 
quark- ant iquark jets (right panel, dash-dotted line). The contribution of the diagram A is 
somewhat bigger than presented in the literature in particular when using CHID matrix 
elements (left panel, short dashed line). The main reason is that in the literature (see e.g. 
Ref. [9}) rather large lower cuts on screening gluon transverse momentum are imposed 
(1-2 GeV). Here we use gluon distributions which allow to decrease the lower cut on the 
gluon transverse momenta in the amplitude down to g™ n ~ 0.4 GeV 2 which is consistent 
with the Tevatron data on exclusive production of \ c data [24]. Having such a low cut in 
the case of dijets production, we rather overestimate the experimental cross section. 

The contribution corresponding to diagram B turns out to be much smaller than that 
for diagram A known from the literature. In addition, it falls much steeper with minimal 
ET,min- In the case of quark-antiquark dijets we present the contribution of uu, dd, ss, cc and 
bb. In the first three cases, we put the quark masses to zero, and in the last two cases we 
take explicit masses known from the phenomenology (1.5 GeV and 4.75 GeV, respectively). 
The sum of all quark-antiquark contributions is shown in the right panel by the dash-dotted 
curve. We conclude that the quark-antiquark jet contribution is smaller by more than two 
orders of magnitude than the digluon one. However, as shown in Ref. the bb contribution 
can be essential e.g. as a background for Higgs searches in exclusive pp scattering. 




FIG. 5: The total cross section as a function of -E^mm- The experimental data points are taken 
from Ref. [l|]. Left panel: digluon contribution for diagam A with our matrex element (solid line) 
and CHID matrix element (short-dashed line), for diagram B (long-dashed line). Right panel: 
quark-antiquark (dash-dotted line) contribution. 



The exceptional dominance of digluon jets over quark-antiquark jets found here offers 
extraordinary conditions for increased glueball production in gluon fragmentation 27] . In 



order to investigate it more one needs to study a contamination of central diffractive com- 
ponents where the proportions of digluonic to quark-antiquark jets are less favourable. 

Similarly to Ref. [9| , we present model uncertainties due to the choice of PDF (FigEJ) and 
due to the choice of the scale (/ip = Hr) in the left panel of Fig. [7] for jet E™ n dependence 
and for dijet mass distribution in the right panel where we show separately uncertainties for 
diagam A and B. We observe that they are much smaller for diagram B. In the latter case the 
lower curve corresponds to fi 2 = Mjj and the upper curve corresponds to fi 2 = (pl ± +pl ± )/2. 
The smaller uncertainty for diagram B can be explained as follows. The gluon propagators 
cause that in a typical situation when two gluons coming from the same proton line are 
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Jet E^ in (GeV) 
FIG. 6: Uncertainties due to the choice of PDF. 

hard and the other two are soft. This is different compared to diagram A where typically 
all gluons are rather soft (or semi-hard). 




JetEf 1 (GeV) M fi (GeV) 



FIG. 7: Uncertainties due to the choice of the scales. Left panel: jet Er distribution for diagram 
A, right panel: invariant mass distributions for both diagrams. Details are explained in the text. 

In Fig. [E] we show pseudorapidity distributions of one of the gluonic jets (left panel) and 
distribution in jet pseudorapidity difference (right panel) for the fixed lower cut on Ej<. The 
distribution for naively calculated diagram B is flatter than that for diagram A. The same is 
true in rapidity difference where the two contributions are almost identical for large rapidity 
differences where, however, the cross section is rather small. The corresponding distributions 
for improved method for calculation of diagram B are quite different than those for naive 
calculation. The contribution of quark-antiquark jets (dash-dotted curve) is negligible. The 
rapidity distributions were not presented by the CDF collaboration in Ref. 

In the left panel of Fig. we show in addition the corresponding distributions in jet 
transverse momentum. The contribution of diagram B becomes negligible at large jet trans- 
verse momenta (or transverse energy). In the right panel we show the distribution in dijet 
invariant mass. At large invariant masses, naive calculations of diagram B give cross sections 
which are similar to the leading contribution from diagram A. However, the contribution of 
diagram B from improved prescription is sizeable only at small invariant masses and does 
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FIG. 8: The distribution in jet pseudorapidity (left panel) and in pseudo-rapidity difference (right 
panel) . 

not have any meaning in the important for Higgs searches large Mjj region. 




P T (GeV) Mjj (GeV) 



FIG. 9: The distribution in jet transverse momentum (left panel) and in dijet invariant mass (right 
panel). 

In Fig. [TU] we show separately contributions for different helicity combinations. Being 

fully consistent with the J z = 0, we see the dominance of the ++ = contributions over 

H — = — h ones. However, our helicities are in the proton-proton center-of-mass system so 
the relation to the J z = rule is only approximate and strictly valid in the high-p^ jets 
limit. 

We show similar distributions at nominal LHC energy yfs = 14 TeV in Fig. [11] and [12j 
The situation is qualitatively similar as for the Tevatron case. Here the distributions in jet 
transverse momentum and dijet invariant mass are somewhat natter. At the LHC energy, 
the contribution of diagram B is much smaller than for diagram A in the whole range of 
kinematical variables considered here. 

Finally, we wish to discuss correlations in rapidity between gluonic jets. In Fig{13] we 
show distribution for diagram A (left panel), naively calculated contribution of diagram B 
(middle panel) and contributions of diagram B calculated as proposed above (right panel). 
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FIG. 10: The distribution in jet transverse momentum (left panel) and in dijet invariant mass 
(right panel) for different dijet helicity states. Only diagram A is included here. 
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FIG. 11: The distribution in jet pseudorapidity (left panel) and in pseudorapidity difference (right 
panel). 

The correlation for the first and second case is similar. The two gluonic jets are only 
weakly correlated in these variables. It is completely different for the contribution of diagram 
B calculated using the ratio of collinear off-diagonal gluon distributions (see Eq. fl4.7p ). We 
observe a strong ridge along the diagonal y 3 = ?/ 4 . This ridge is a consequence of the ratio 
defined in Eq.fl4~Tl). 

In Ref . 0| we have studied in detail irreducible exclusive bb background to exclusive Higgs 



boson production. The gluonic jets can be misidentified as 6-quark jets |28j. If both gluonic 
jets are misidentified then such a misidentified event can contribute to a background to 
exclusive Higgs boson production. In Fig[JJ]we illustrate the situation. We show both the 
Higgs signal (hatched area) including experimental resolution 29, 3(| as well as diffractive 
bb continuum, QED bb continuum as well as formally reducible digluon contribution. In the 
calculation we have assumed that jet misidentification probability is 1.3% (which corresponds 
to the ATLAS misidentification factor with a b-tagging efficiency of 60% j28l), i.e. we have 
multiplied the dijet cross section by a quite small number 0.013 2 . The obtained contribution 
is even larger than the bb one and overlays the Standard Model Higgs signal. In the case of 
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FIG. 12: The distribution in jet transverse momentum (left panel) and in dijet invariant mass 
(right panel). 



diagram A 



diagram B: naive prescription 



diagram B: ERBL >. - 0.4 






FIG. 13: da/dysdy^ for diagram A (left panel), for diagram B, naive prescription (middle panel) 
and for diagram B calculated in the way proposed in this paper (right panel). 



exclusive production of Higgs boson beyond the Standard Model the situation can be better 
(see e.g. Refs. 0, M, H^)- 



VI. CONCLUSIONS 



In the present paper we have discussed the exclusive production of dijets (both digluon 
and quark-antiquark ones). We have included the contribution, previously known from the 
literature, when both gluons are emitted from the same t-channel gluon line, as well as 
new contributions when the gluons are emitted from different t-channel gluon lines (see 
Fig. [T]). We have presented corresponding formulae with simple prescriptions for the unin- 
tegrated gluon distributions relevant for these two cases. For both contributions, we made 
predictions for various differential distributions at Tevatron (1.96 TeV) and LHC (14 TeV) 
energies, including an analysis of different gluon polarisation contributions and theoretical 
uncertainties. 

The diagram B contribution turned out to be much smaller than the one (diagram A) 
known from the literature. They become comparable only for large rapidity differences of 
both jets where the cross section is rather small or when the jet transverse momenta are 
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: p p -> p p b b \fs = 1 4 TeV : 




M bB (GeV) 

FIG. 14: Invariant mass distribution of the bb system. Shown are contributions from diffractive 
Higgs boson (shaded area), bb continuum (solid line), 77 continuum (dash-dotted line) and diffrac- 
tive digluon contribution (dashed line) multiplied by an ATLAS misidentification factor squared. 

small. The latter case can be, therefore, very important for the central diffractive production 
of pions. This will be discussed elsewhere. 

We have found that jets corresponding to the mechanism of diagram B are strongly 
correlated in their rapidities. This is a consequence of a specific behavior of off-diagonal 
gluon distributions in the ERBL region where \x\ < |£| <C 1. 

We have compared our results with the CDF collaboration data. Our cross sections are 
somewhat larger than those obtained in the literature. Compared to those calculations, we 
have performed integration starting from smaller lower limit for screening gluon transverse 
momenta and with the choice of scale in the Sudakov form factor as advocated by Coughlin 
and Forshaw [3|. Our observation may mean e.g. that the gap survival probabilities are 
smaller than usually assumed for this process. Clearly, further work on this issue is required. 

We have discussed also the dijet reducible background to the central exclusive Higgs boson 
production. In the framework of the same model, we have calculated the contribution of the 
quark-antiquark jets CEP. We have found that this contribution is much smaller than that for 
the gluon-gluon dijets in the whole phase space. However, the quark contribution, especially 
the bb one, is very important, as it constitutes an irreducible background for exclusive 
production of the Higgs boson. On the other hand, the gluonic jets can be misidentified as 
the bb jets and in this sense they also contribute to a background for exclusive production 
of the Higgs boson. 

If both gluonic jets are misidentified as b or anti-6 quark jets, then this leads to an extra 
background to exclusive Higgs boson production when the Higgs boson is observed in the 
bb decay channel. This extra contribution can be even more important than the irreducible 
bb contribution. When the realistic ATLAS misidentification factor is included one obtains 
the total background which significantly exceeds the Higgs signal. This observation suggests 
that the experimental observation of the exclusive Standard Model Higgs production may be 
very challenging. The situation may be better for beyond the Standard Model Higgs boson 
production though, but corresponding detailed analysis including a Monte-Carlo simulation 
of backgrounds still needs to be done. 
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